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Abstract 

We construct a chiral gauge theory to describe fractionalization of fermions in graphene. 
Thereby we extend a recently proposed model, which relies on vortex formation. Our chiral 
gauge fields provide dynamics for the vortices and also couple to the fermions. 

1 Introduction 

In some condensed matter systems the electron excitation spectrum near the Fermi surface can 
be described by a Dirac-type matrix equation. This equation does not arise from relativistic con- 
siderations, but rather by linearizing the energy dispersion near (a finite number of) Dirac points 
(intersections of the energy dispersion with the Fermi level). Such systems can exhibit fermion 
fractionalization if the Dirac equation possesses isolated bound states in the gap between negative- 
energy (valence band) states and positive-energy (conduction band) states. 

A familiar example is 1-dimcnsional polyacetylene [Tj. As is generally the case in one dimension, 
there are two Dirac points for polyacetylene. Therefore a 2-component Dirac equation governs 
electron motion near each Dirac point. A distortion of the underlying lattice (Peierls' instability) 
perturbes the electron motion in a way that couples the two Dirac points and opens a gap in the 
energy spectrum. In the Dirac equation description this is achieved by coupling the Dirac field 
4" to a scalar field <p, which is a measure of the lattice distortion, tp enjoys a symmetry with 
two ground states in which it takes homogenous values ± (fQ. This coupling to ipo leads to a Dirac 
mass tx \ipo\. But tp can also take a position-dependent kink profile (soliton) tp s that interpolates 
between the two vacua ±<po- This "twisting" of the mass parameter describes a defect in the lattice 
distortion. The Dirac equation with the kink profile tp s replacing the homogenous ±<po possesses 
a single mid-gap (zero-energy) eigenstate. This gives rise to fractional fermion number for the 
electrons: 1/2 per spin degree of freedom [T]. 

Recently a similar story has been told by C.-Y. Hou, C. Chamon and C. Mudry [2] (HCM) 
about (monolayer) graphene. This is a 2-dimensional honeycomb array of carbon atoms forming a 
hexagonal lattice, which may be viewed as a superposition of two triangular sublattices, A and B. 
The generators of lattice A are ai and a2 . The three vectors Si connect any site from lattice A to 
its nearest neighbor sites belonging to B. They are 



Bl = (0,-1)4 Ba = * -U s 3 = ~.o K (!) 
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Figure 1: Graphene hexagonal lattice constructed as a superposition of two triangular lattices A 
and B, with basis vectors ai for lattice A and vectors Si connecting A to B. 



where t is the lattice spacing; see Figure 1. 

When no lattice distortion is considered, the tight-binding Hamiltonian, with uniform hopping 
constant t, is taken as 

F « = -*E E (aHv)b(r + Si )+b\v + Si )a(r)\, (2) 

reA i=l,2,3 ^ ' 

where the fermion operators a and b act on sublattices A and B. Hq is diagonal in momentum 
space. 

a(k)l_ Ee _, k . r fa(rjl (3) 



b(k) ^ \b(r 



reA 



H ° = Y1 ( $ ( k ) at ( k ) b w + ( k ) 6t ( k ) a ( k )) 



(4a) 



$ (k) = -t e *' s< ( 4b ) 

i=l,2,3 

The single particle energy spectrum E(k) = ± | $(k) | contains two zero-energy Dirac points at 

k = K± = ± G^l? ) ; * <K±) ^ a 

This parallels the 1-dimensional case, but is rare in two dimensions J3J. 

Ho is linearized around the two Dirac points, k = K± + p, and is supplemented by a term 
arising from a (Kekule) distortion of the lattice, which couples the two Fermi points. 



H = E ( ^+ (p) a +(p) b + (p) + ^+ (p) b + (p) a + (p) + ^- (p) a - (p) 6 - (p) + ^* (p) &+ - (p) a - (p) 
p ^ 

+ ]T f A a^ (p) 6- (p) + A* bl (p) a+ (p) + A* a + _ (p) b + (p) + A 6+. (p) a_ (p)) (5) 
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<j>± is the linearization of $ : 4>±(p) = ±vp (p x ± ip y ),VF = 3t£/2, (hence forth vf is set to unity) 
and a±, b± are fermion operators near the Dirac points: a±(p) = a (K± + p), b±(p) = b(K± + p). 
In the second sum Ao is a homogenous, but complex order parameter, which effects the coupling 
between the Fermi points K± and leads to a mass gap in the single-particle energy dispersion: 
e(p) = ±v/p 2 +l A | 2 . 

To find zero modes for this system HCM promote the mixing parameter Ao to an inhomogenous 
complex function A with a vortex profile. To this end the Hamiltonian (5) is presented in coordinate 
space as 

H = J d 2 r^* (r) A'*(r), (6) 
where "J (r) is a 4-component "spinor" 

V 4t 



with V± (r) = / d 2 pe- l P r a± (p) 
i>± ( r ) = / <Ppe- ip - T b± (p) 



and K is the 4x4 matrix 



K 





-2id z * 
A*(r) 




-2id z 



A*(r) 



A(r) 



2id z * 





A(r) 
2id z 




with -2id z = \ (d x -id v ). 

HCM take A(r) in the n- vortex form: A(r) e 1 



where n is an integer, A(r) vanishes as 



(7) 



■ I' 



for small r, and approaches the mass-generating value Ao at large r. They then establish the 
occurrence of |n| zero modes, i.e. solutions to K^> = 0, on lattice A (B) for negative (positive) n, 
and they construct explicitly the solution for n = — 1. Therefore fermion number oc J d r^*^f is 
fractionalizcd. 



2 Chiral Gauge Theory for Graphene 

In this paper we elaborate the HCM model, and address the following two topics. HCM leave un- 
specified the dynamics that gives rise to the complex vortex profile. To remedy this, we first propose 
identifying the HCM vortex with the Nielsen, Olesen/Landau, Ginsburg, Abrikosov (NO/LGA) vor- 
tex, which is described by a charged scalar field, as in the HCM model. However an U(l) gauge 
field is also involved in creating a finite-energy NO /LGA vortex, but no gauge field is present in the 
HCM model. Therefore, second, we propose introducing a relevant gauge potential and coupling 
it to the Dirac fermions in a chiral manner. This expands the symmetry of the interaction to the 
kinetic term, and renders the theory invariant against local, chiral E/(l) gauge transformations, 
which also act on the scalar and Dirac fields. Additionally there remains a global fermion number 
U(l) symmetry, and its charge is fractionalizcd. 

To present our extension of the HCM model, we shall use Dirac matrix notation. We begin by 
rewriting K in (7) in terms of Dirac matrices whose forms we take as follows. 

Note: we use 4x4 Dirac matrices, even though the minimal Dirac algebra in (2+l)-dimcnsions 
requires only 2x2 matrices. But we have four degrees of freedom: two each in lattices A and B. 
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Since there are two spatial dimensions, we use only the first two a matrices : a l ,i = 1,2 or x, y; 
while the role of the third one, a 3 , will emerge later. Also unlike with minimal 2x2 Dirac matrices, 
here we can construct the chiral 75 matrix, as the Hcrmitian quantity 



The 7 matrices read 



1 I I 

75 = -101 a a = ( Q _ J ) , 7 5 = /. 



7 = /3a=r° f), 7° = /?, 75=*7°7 1 7 2 7 3 



<X 



Note that the a 1 (i = 1,2,3) anti-commute among themselves and with (3; 75 commutes with the 
a 1 and anti-commutes with (3. 

With these matrices K in (7) may be presented as 

A'* = (a ■ p + g/3 [</ - itp 1 7,5]) *. (8) 

(Henceforth all vectorial quantities are 2-dimensional.) Here p is the operator — iV; we have 
renamed A as g tp, where g is a coupling strength and tp is a complex scalar field, with real and 
imaginary parts: tp = tp r +itp l . Note that when (p is decomposed into modulus and phase: tp =\ 
tp I e lx , the interaction part of © may be presented as g \ ip \ ^* (3 e~ %l5X ^ . This makes it clear 
that this interaction is invariant against a local chiral gauge transformation. 

tp -> e 2tuJ tp^> x^X + ^u 

[When (ys = 920 is constant, its constant phase may be removed from K by the above transformation, 
leaving a conventional Dirac mass term (gap) oc \tpo\. ] 

In order that the kinetic portion of ([8]) be invariant against the local gauge transformation ([9]), 
we introduce coupling to a gauge potential A, which transforms as 

A^A + Vcj. (10) 

Thus our final Dirac Hamiltonian density reads 

Ka * = ** a ■ [p - 75 A] ^ + g^*f3 [<p r - 175 99*]* 

= 7 ■ (P - A) + 7 ■ (p + A)*- + g ip §+ + g^* . (11) 

We have introduced the Dirac adjoint ^ = 7 , and the chiral components *S?± = \ (1 ± 75) 'if, 
whose gauge transformation law is 

# ± e ±iw *±, * ± ^* ± e T ^. (12) 

The Bosc fields tp and A arc determined by the familiar NO/LGA equations. 

D • D tp = tpV (tp* tp) 

T>tp= (W - i2A) ip (13) 

^B ii d j B = 3 i OBR (14) 

B ee e « i), A j 

Jbose = 4/m tp* Dtp 
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Hero e is a further coupling constant and V is chosen so that at minimum V' = 0, <p* tp = ip^, 
which gives rise to the Dirac mass for $ cx \ipo\. Furthermore, vortex solutions to these equations 
also exist. Their form is 

ip v {r) = <p{r)j nB , 

A;(r) = -n f '^a(r), (15) 

where n is an integer; <p(r) vanishes with r as r' n ' and tends to ip$ at infinity; a{r) vanishes at the 
origin so that A t , is regular there, and a(r) tends to 1/2 at large r. All this ensures finiteness of 
the vortex energy J d 2 r B 2 + | D tp | 2 + V ((p*tp)). 

Finally, note that our system possesses a global fermion number symmetry, with just the Fermi 
fields transforming with a constant phase: $ — > e lX \&. Consequently the theory possesses a local 
chiral U(l) symmetry and a global U(l) fermion number symmetry Because the theory resides in 
(2+1) dimensions, no chiral anomalies interfere with our chiral gauge symmetry. 

Indeed, in spite of the presence of the 75 matrix, the theory is parity (P) and charge- conjugation 
(C) invariant and the axial vector current is C-odd. This is because in (2+1) dimension, with 4- 
component Dirac fields, the relevant transformations read 

P : cp (t,x,y) -> ip (t,-x,y) 

A°>y {t,x,y)^A°>v {t,-x,y) 

A x (t, x, y) -> -A x (t, -x, y) 

^{t 7 x,y)^i 1 3 1 1 ^ (t,-x,y) (16) 

C : (p — > ip* 
A -> A 

- 7« *j (17) 
It follows that time reversal symmetry holds also. 



3 Modified Dirac Equation 

With the additional gauge potential A, our Dirac eigenvalue problem differs from HCM. According 
to (fTTj) we have 

(a ■ (p - 75 A) + gp[<p r - * 7s ^*])* = £? * . (18) 

Observe that a 3 , which we rename R, anti-commutes with the matrix structure on the left side of 
(fl"8]) . Therefore if is an eigenfunction with eigenvalue E, R^e belongs to eigenvalue —E, and 
zero modes can be chosen as eigenstates of R. This is a consequence of the "sublattice symmetry" 
identified by HCM. 

Next we show that the gauge interaction in (fT5|) does not affect the zero modes found by HCM 
without A. To this end, we adopt the Coulomb gauge and present A as A 1 = e y dj A. Also it is 
true that a 4 75 = —ie %3 a?R (i = 1,2). Consequently the kinetic term in fTS"]) can be written as 
e~ AR a ■ p e~ AR , and (JTSJ) becomes 

(a.p+$/% r -*W]) (e- AR *)=E (e AR y). (19) 

Thus e~ AR VP satisfies the HCM equation at E — 0. Comparison with (fT5"|) shows that A'(r) = 
—n a (r) / r, so that at infinity A (r) tends to — ^ In r, and the zero modes with the gauge interaction 



■5 



differ at large r from the HCM modes by factors y-^"/ 2 ) . This does not affect nomalizability because 
the zero modes found by HCM are exponentially damped by the interaction with ip v . Finally, since 

/ \ 



the n = — 1 HCM mode as well as ours has the form, 



(' 

V o / 



where v is an exponentially 



damped function, we see that indeed is an R eigenstate, with eigenvalue —1. Fermion number 
fractionalization in the gauge theory is now established by the same reasoning as in HCM. 

4 Energy Relations 

The total energy functional for all our fields is 

£total = J d2r {j^ B 2 + I | 2 +V(y>»j + j d 2 rV*K A «f. (20) 

Varying this with respect to \&* produces our Dirac equation (|18|) at zero eigenvalue. Varying with 
respect to the Bosc fields ip* and A derives (fT3|) and (fl4|) . but with a back reaction from the Dirac 
fields. 

T>-r>v = vV' + | * (1 + 75) * (21) 

^'8/5=^ + $^$ (22) 

However, with a zero mode that is an eigenstate of R, the back reaction Dirac bilinears vanish. 
Thus our Dirac zero mode, together with the scalar field/gauge field NO/LGA vortex, is a self 
consistent solution of the coupled system. 

Chiral gauge theories have previously entered physics, but in even-dimensional space-time, where 
the chiral anomaly influences the structure and physical utility of these models j4j . In the present 
work, we have a chiral gauge theory in odd-dimensional space-time, whose structure is mathemat- 
ically very elegant owing to its self-consistent solutions. It remains to be determined whether a 
microscopic description for graphene can lead to the chiral gauge field that enters our theory. 
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